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O ' Abstract 

o 

CN ■ We consider the Schrodinger operator with a periodic potential on a quasi ID contin- 

uous periodic model of armchair nanotubes in M 3 in a uniform magnetic field (with 
amplitude B 6 1), which is parallel to the axis of the nanotube. The spectrum of 



< 



this operator consists of an absolutely continuous part (spectral bands separated by 
gaps) plus an infinite number of eigenvalues with infinite multiplicity. We describe all 
eigenfunctions with the same eigenvalue including compactly supported. We describe 
the spectrum as a function of B. For some specific potentials we prove an existence of 
gaps independent on the magnetic field. If B ^ 0, then there exists an infinite number 
of gaps G n with the length \G n \ — > oo as n — > oo, and we determine the asymptotics 
of the gaps at high energy for fixed B. Moreover, we determine the asymptotics of the 
gaps G n as B — > for fixed n. 



1 Introduction and main results 



o 
o 



We consider the Schrodinger operator M'b = (— iV — £/) + V q with a periodic potential V q 
on the armchair nanotube T N C M 3 , iV ^ 1 in a uniform magnetic field SB = B(0, 0, 1) G M 3 , 
Bel. The corresponding vector potential is given by srf{x) = ^[3§,x] = |(— x 2 , Xi, 0), x = 
(x\, x 2 , x 3 ) G M 3 . Our model nanotube T N is a union of edges r w , i.e., 

O ' 

r w = U we2 r w , u = (n,j,k)eZ = ZxN 6 xZ N , R m = {l,2,...,m}, Z N = Z/(NZ), 

see Fig. (TJEJ Each edge = {x = r^ + te^, t G [0, 1]} is oriented by the vector e w G IR 3 and 
has starting point G M 3 . We have the coordinate x = r w + te w and the local coordinate 
t G [0, 1] (length preserving). We define r u , e^, u> = (n,j, k) G Z by 

r n ,i,fc = 2n/ie 3 + -R(c 2fc , s 2 &,0), r ni4jfc = 2nhe 3 + i?(cos/5 fc , sin/5 fc , 0), /3 fc = 2(3 + 2 fc, 
r n ,2,fe = r„, 5 ,it = (2n + l)/ie3 + i^(cosCfc,sinCfe,0), ( k = (3 - a + 4> 2kl </>fc = — , 
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Figure 1: Armchair graph for N = 10 and for N — 1. 
r„, 3 ,fc = r nAk = (2n + l)he 3 + R(c 2 k+i, s 2 k+i, 0), 

1 l Jcos^ + I 

e„, 6 ,jfc = r n+1(4 ,jfe - r n ,6,jfc, sin/3 = —, sina = — , R= — , (1.1) 

R 2R sm0i 

0<ct</3,a + /3 = fand/i = V2 + i2i-R 2 - 2i? 2 , i? p = ^/(p.R) 2 - l,p = 1,2. 

Each edge is the segment and connect the vertices, which lie on the cylinder C = 
{(xi,x 2 ,Xs) G M 3 : x\ + x\ = R 2 }. Note that each edge (without the endpoints) lies 
inside the cylinder C. The points i"o,i,fe are vertices of the regular N-gon = 1,2,3,4. 

■^3 (^2 ) arises from &\ (^4 ) by the following motion: rotate around the axis of the 
cylinder C by the angle and translate by he 3 . Repeating this procedure we obtain T N . 

For each function y on T N we define a function y u = y\r u ,^> G Z. We identify each 
function y^ on with a function on [0, 1] by using the local coordinate t G [0, 1]. Define 
the Hilbert space L 2 (T N ) = ®^zL 2 {T u ). Our operator M'b on T acts in the Hilbert space 
L 2 (T N ) = © w L 2 (r w ) and is given by 

= -5£/o,(*) + q(t)f u (t), d 0J = --ia UJl a w {t) = (^(r w + tej, ej, (1.2) 

see jSDD] . where = 9 / w , g G L 2 (0, 1) and © w /^' G L 2 ^) satisfies 

Magnetic Kirchhoff Boundary Conditions: y is continuous on T N and 

du 2 yu, 2 (0) - d^y^il) + d^y^O) = 0, d^y^O) - d^y^il) + d^y^O) = 0, 
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Figure 2: A piece of a nanotube r . The fundamental domain is marked by a bold line. 

-d U3 y U3 (l)+d u)4 y tJA (0)-dy n - 1A k{'i-) = 0, -d Ui y UA (l) + dy n ^ k+1 (0) - dy n - li5>k +i(l) = 0, 

(1-3) 

for all ujj = (n,j, k), (n,j, £;) G Z x N 6 x % N . 

Condition (pQ) means that the sum of derivatives of y at each vertex of T N equals and the 
orientation of edges gives the sign ±. Such models were introduced by Pauling |Pa] in 1936 to 
simulate aromatic molecules. They were described in more detail by Ruedenberg and Scherr 
[RSJ in 1953. Further progress had been made toward periodic systems by Coulson in [Co] 
where a network model of graphite layer was worked out. A network model of a crystal with 
non-trivial potentials along bonds was studied by Montroll [Moj . Carbon and boron nitride 
nanostructures, in particular nanotubes, graphene layers have been modelled by quantum 
networks with a honeycomb lattice structure in [ALMj . |ALM1] . As it is shown in [HaJ, 
Ch.7.3.3, the structure of carbon-boron-nitride single wall nanotubes may be complicated. 
The edge potential q is not even for such tubes. We consider one of nanotube models, where 
the potential is not even. The other models can be considered by the similar way. 

The standard arguments (see [KLj ) yield that ,Wb is self-adjoint. For simplicity we shall 
denote r^i C T 1 by r a , for a = (n,j) G Z\ = Z x Ng. Thus T 1 = U a eZi^a, see Fig CD In 
Theorem 1 1.1 1 we will show that M'b is unitarily equivalent to H a = (BiH^, where the operator 
El acts in the Hilbert space L 2 (T l ) and is given by (H%f) a = -/„+?/«, (f a )ceZi, UDaeZi e 
L 2 (r 1 ), and the components f a satisfy: 

e mi /n,l(l) = /n )2 (0) = /n,5(0), e M2 /n,2(l) = /n, 3 (0) = / n , 6 (0), 

e Wl /„, 3 (l) = /n, 4 (0) = e iai / n _ li6 (l), e ia *s k f nA {l) = f n>1 (0) = / n _ 1)5 (l), (1.4) 
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e mi /; 3 (l) - f'n,M + e i0l /;- M (l) = 0, e ia *s k f nA {l) - /^(O) + e"^^_ li5 (l) = 0, (1.5) 
where 



2iri 



^_ B{R 2 — Ri) BR 2 B(2R 2 — Ri) , . 

e", a! = , a 2 = , a = ax + a 2 = , (1.6) 

R is given by (11 .11) . and R v = \J (pR) 2 — l,p = 1, 2. For the operator H£ we define funda- 
mental solutions tfjj. 1 ^ = (^jce^, yjj^ = (^) ae ^ u z/ = 1,2 which satisfy 

-Ay + V q y = \y, on T 1 , b. c. (HI, (USD 

e M (l ) A) = $ / M (l,A)=J 2 , e' M (l,A) = $ M (l,A)=0, (1.7) 

where /„, n ^ 2 is the n x n identity matrix and 



\qW ,q(2) \ /.-(i) .-(2) 1 

i fc,n,6 fc,n,6/ \r 6 rfc,n,6. 



We define the monodromy matrix by 



m ' (a)= k; c; (1 ' a) - (L9) 

This matrix is determined by the fundamental solutions on a fundamental cell Tq = U^^Tqj. 
In order to analyse the operator H a we use the methods based on a research of the mon- 
odromy matrix M k (see [EE], [KL1] . |BBKL] . |BBK1] . and also [BBK] . [UK] . [E])- 

Let cr£) = {fi n ,n ^ 1} be the spectrum of the problem — y" + qy = Xy,y(0) = y{l) = 
(the Dirichlet spectrum). We formulate our first result. 

Theorem 1.1. (i) The operator Jrf? B is unitarily equivalent to H a = ®^H%. 
(ii) Let (a, fc) e 1R x Zjv- Then for any A G C \ od there exist unique fundamental solutions 
^ = ($ki)aez 1 ,<Pk = [wfyaez^v = 1,2. Moreover, each of the functions $^ a (x,\), 
Vk(J\ x -> A), x G T 1 is analytic in A G C \ a~£> and the monodromy matrix -Mjt(A), A ctd has 
eigenvalues of the form t^\, and satisfies 

detA^ = l, M k (X) T JM k (X) = J, where J=(°. l) , j = (? lV (1.10) 



-j oy ' J 

Furthermore, the matrix-valued function TtM. k Tlr x is entire, where TZ = I 2 © <^i-^2- 

For all S G R the monodromy matrix A4 k has a simple pole at each point A G o~z>, which is 
an eigenvalue of H% (see Theorem I2.2p and has no any other singularities. Moreover, these 
poles are independent on the magnetic field. This result is in contrast with the corresponding 
result for the zigzag tube (see |KL] , |KLlj ). For the zigzag tube there exist some singular 
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values of the magnetic field, where the monodromy matrix is not well defined. Moreover, the 
spectrum is pure point for such singular values and the spectral bands shrink to the points 
as the magnetic field approach to the singular value. 

For the equation —y" + q{t)y = Xy on the real line we introduce the fundamental solutions 
0(t,A) and (p(t,X),(t,X) 6lxC satisfying 0(0, A) = ^'(0,A) = 1,0'(O,A) = <p(0, A) = 0. 
We define the functions F, F_ by 

F=^p, F_ = ^-^, (1.11) 

where (fx = ip(l, -),0i = 0(1, •), = <f'(l, ■),$[ = 0'(1, •)■ We formulate the results about 
the Lyapunov function and the spectrum of H k . 

Theorem 1.2. Let (a,k) Glx 7L N and let r k \ , r k \ be eigenvalues of M. k . Then 

(i) The Lyapunov functions F kjV = \{jk,v + T kl)i v = 1,2 are branches of F k = + J~P~k on 

the two sheeted Riemann surface fR k defined by J~pk and satisfy: 



9i?2 _ p2 _ 2 

f m = ^ - (-irv^, & = — if — ^ pk = ( 9f2 - s I)4 + sIf 2 _, (i.i2) 

det(M fc - r/ 4 ) = (r 2 - 2F M r + l)(r 2 - 2F fej2 r + 1), (1.13) 

where 

Tck irk 
c fc = cos(— + a), Sfc = sin(— + a), (1.14) 

a zs given by $1.6}) . 

(ii) If -Ffc(A) G (—1, 1) /or some A G K and A is not a branch point of F k (X), then F k (X) ^ 0. 
(Hi) The following identities hold: 

a{H a k ) = a pp (H a k ) U a ac (H a k ), a pp {H a k ) = a D , 

a ac {H a k ) = {X G R : F kjV (X) G [-1, 1] for some u G N 2 }, (1.15) 

a{H k a ) = a(H a N _ k ), a{H?*) = o{H a k+l ), a(H^) = a(H a k ) (1.16) 

(each point of the spectrum is counted with multiplicities). Each point Ao G a pp (H k ) has an 
infinite multiplicity. 

Remark. 1) We take the branch of ^fpi such that if Pk(X) > 0, A G R, then a/ Pk(X) > 
and = £ k + ^/p^ > F k 2 = C,k ~ \fp~k- Then ^/p^ is not an entire function even for c k = 
or s k = 0. Using other branches we could define a new function ^/p^ which is entire for 
c k = or s k = 0, but this choice is not convenient for labeling of endpoints of gaps (see 
below Theorems 1 1 . 34TT761) . 

2) The statement (ii) is similar to the standard result of Lyapunov for the Hamilton systems. 
The similar results hold for the zigzag tube in magnetic field (see |KL] . |KL"T] ) . 

3) Periodicity of the spectrum with respect to B holds for the zigzag nanotube (see [KLlj ). 
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rpk,+ rpk- 771O,- 
-^1,2 -^1,3 -^1,3 



7-,/c ,- J-,*,- 
-^2,2 -^2,2 



-^2,3 



77.U- T7IU- R U,-h TTlU.-f 

'2,0 ^2,1 ^1,1 ^1,1 -^2,1 

Figure 3: Graph of the function F k (X) and the spectrum of H k 



v 2,3 



4) Below we will sometimes write c k (a),F k}U (X,a), ... instead of c k , F k>1/ (X), when several 
magnetic fields are being dealth with. In order to prove (I1.16P we use the identities 

7T 7T 

sjfc(-a) = s N _ k (a), c fc (-a) = c N „ k (a), s k (a + —) = s k+1 (a), c k (a + — ) = c fe +i(a). (1.17) 

5) We say that {A } C a pp (H^) is a flat band. The corresponding eigenfunctions of H% have 
compact supports (shown by FigHJ) and are described in Theorem 12.21 

We define the entire functions 



&et{M k =F h) = 4(F M t l)(F k , 2 =F 1). (1.18) 
The zeros X*'^ n ,n ^ 0, v = 1,2, of the function are the periodic eigenvalues. The zeros 



A 



k,± 
u,2n- 



. x , (z/, n) 6 N2 x N, of Z) fe are the antiperiodic eigenvalues. Let 



A 



2,0 



\K,T ^ \k,- ^ \K,~ ^ \K,-T ^ \K,~ ^ \K- \K,^- ^ >«,t ^ 

^ A l,0 ^ A l,2 A 2,2 ^ A 2,2 A l,2 ^ A l,4 A 2,4 ^ A 2,4 A l,4 ^ 



and 



A 2,l ^! A l,l ^ A l,l ^ A 2,l ^ A 2,3 ^ A l,3 ^ A l,3 ^ A 2,3 ^ ••• 

counted with multiplicities. Such labeling was introduced in jBBKl] for the case a = 0, and 
is associated with the Lyapunov functions F ki i,F kj 2 (see Figj3j). 

A zero of p k , k e is called a resonance of H k . There exist real and non-real resonances 
(see example from |BBKLj ). Roughly speaking the simple real resonances create gaps. Note 
that in the case of zigzag nanotubes all resonances are real (see |KLj . |KLlj ). 

We define the functions 



u k (X) = \F„(X)\ - si, v h {\) = \F.{\)\ 



x e 



k e Zw. 



;i.i9) 
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Theorem 1.3. Let (a, k) G R x Zjy. 27ien the identity 
holds, where each spectral band satisfies: 



^, P -i = •Vp-ij E 2, P = A 2iP , E 1<p = \ * ' p ± p = 2n - 1, (1.21, 



A^(a)=A° : J(0), JSfe ± (a) = J Eg : (0), x n (a) = x n (0), (1.22) 
u>/iere 



fc,n 



min{A G x n : p fc (A) = 0}, r^ n = max{A G x n : p fc (A) = 0}, x n = (X°{ p , A?'* ). 

'(1.23) 



Remark. 1) For some specific potentials and some n ^ 1 the intervals x n are antiperiodic 
gaps (see below Proposition 13.51 and Theorem 11.61 (i)). These gaps are independent on the 
magnetic field. 

2) We describe a ac (H%) in details in Theorems 13.41 [3T6l 

Recall the needed properties of the Hill operator Hy = —y" + q(t)y on the real line 
with a periodic potential q(t + 1) = q(t),t G R. The spectrum of H is purely absolutely 
continuous and consists of intervals [A^L^A"],^, ^ 1. These intervals are separated by the 
gaps 7 n = (A~,A+) of length |7 n | ^ 0. If a gap 7„ is degenerate, i.e. \j n \ = 0, then the 
corresponding segments merge. The sequence Aq < A|f ^ A^ < ... is the spectrum of the 
equation — y" + qy = Xy with 2-periodic boundary conditions, that is y(t + 2) = y(t),t G R. 
Here equality A~ = A+ means that A^ is an eigenvalue of multiplicity 2. It is well-known 
that fi n G [A~,A+],n ^ 1. The function F has only simple zeros r] n ,n ^ 1, which satisfy 
r)i < T]2 < ••■ We describe the spectral gaps of H%. 

Theorem 1.4. Let (a, k) G R x Zjy. T/ien 

ft) cr ac (H%) = R \ Un>oGJ )n , w/iere t/ie #aps G£ n satis/?/: 

70 C Gfc i0 = (-OO, £2,0"), 7n C G^ 4n = (£? 2 ,2nJ ^2,2n)> ^fc,4n-2 = C^l.W-lJ ^l,2n-l) C X «. 
^fc,4n-3 = (-^2,2n-l> ^l,2n-2)> ^t,4n-l = (-^l,2n> ^2,2n-l) 5 £ [-^l,2n-l J -^l,2n-l] • (1-24) 

If sl(a) < sf(ai) for some (ax,£) G R x Z^, t/ien 

^fc,4n C ^2,4n> ^,2n-l C ^fc,2n-l °^ J1 ^ 1. (1.25) 

Moreover, for all n ^ n and some n ^ 1 

G^4n-2 = ^ S ^ ^ 2 , Q 2 n-1 = 0, («, *0 ^ (0, 0). (1.26) 

^ 0, i/(a,A;) ^ (0,0) 



(m,) J/nfc = |F_(A)| — c| ^ on x n /or some n ^ 1, i/ien G£ 



4n-2 



x„ 



We formulate our main result about the spectrum of H a . 
Theorem 1.5. Let a G [0, ^]. Then the following identities hold true: 

a(H a ) = a pp {H a ) U a ac (H a ), a pp (H a ) = a D , a(H a ) = a{H- a ) = a(H a+ %) (1.27) 

(each point of the spectrum is counted with multiplicities). 

(i) a ac (H a ) = R \ U n ^oG a n , where the gaps G a n = ^ k & N G a Kn = (E~, E+) satisfy: 

In C Gl n , G\ n _ 2 C K n , r] n G [E^ n _ 2 , E An _ 2 ], (1.28) 

Gt n GG2, G^CG^, o^a< ai ^^. (1.29) 
Moreover, for all n ^ uq and some n$ ^ 1 

Gtn-2=\ h ] (h > ^ = 0, if (a,N)^ (0,1). (1.30) 

(ii) The spectrum a(H a ), a G (0, ^), has an infinite number of gaps G 2n ^ and \G 2n \ — > oo 
as n — > oo. Moreover, the following asymptotics hold true: 

/ \ 2 

Etn-2m(°>) = \n(n--)±e m ) + q + O(n~ 1 ), n -> oo, m = 0,l, (1.31) 

where go = J q(t)dt,0o = arccos(|y / 5 + 4| cosaj), 6*i = arcsin(|sina). Furthemore, if N ^ 
2 ; i/ien a// ot/ier gaps wi/i sufficiently large energy are empty. 
(Hi) If \F-\ ^ cos 2 a on x n /or some n ^ 1, then G2 n _ 2 — K n- 

Remark. 1) Identities ( 11.27ft show that for analysis of operator H a , a G M. it is sufficiently 
to study a G [0, 

2) Using the statement (iii) we prove in Proposition 13.51 that each G\ n _ 2 = G\ n _ 2 ^ 0, 
(a,n) gMx N no for any no ^ 1 and some non-even potential (which depends on no)- Note 
that if the potential is even, then all the gaps G a An _ 2 depend on the magnetic field (see below 

QEE}). 

Recall the definition of the effective masses for the Hill operator. The identity F(X) = 
cosfc(A), A G C + defines the quasimomentum k(X), which is an analytic function on A G C + . 
With each endpoint of the gap 7„ ^ 0, we associate the effective mass by (we take 

some branches k such that k(X) — > im as A — > A*) 

A = A^ + (fc(A 2 ) M ^ )2 (l + o(l)) as A-A£. (1.32) 

(see |KK] ) . We formulate results for even potentials. Let L^ en (0,l) = {q G L 2 (0, 1) : 
g(l-t) = ?(t),te(0,l)}. 
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Theorem 1.6. Let q G L 2 even (0, 1), (a,n) G [0, ^] x N. Then 

(i) If c k = 0, then p k = and G kAn _ 2 = x n . If c k ^ 0, then all resonances are real and 
G kAn~2 = ( r k,m r tn)- For each fceZ ^ the 9 a P s G a k4n _ 2 satisfy: 

GIm - 2 = * ? 0, ^if {a {a,k) j (0 0) ' Ghn ~ l = ^ (L33) 

If s 2 k (a) < sj(ai) for some (ai,£) G K x 1, N , then G kM __ 2 C G a ^ An _ 2 . 

(ii) The gaps of H a satisfy 



71 



G a in - 2 CGZ_ 2 , 0^a<a 1 ^— , (1.34) 

G ln = In, G\ n _ 2 = r </fl = "° , ^ = 0. (1.35) 

7^ 0, «/ a ^ 



Moreover, if a & (0, t/ien G^n 7^ 0? an< ^ 



£±Ja) = A± + -^ + 0(a 4 ), 7n ^ 0, (1-36) 

n 

Et n (a) = \+± -^t= +0(a 2 ), 7„ = 0, (1.37) 



as a — ► 0. 



Schrodinger operator on graphs is a subject of intensive studies in the last years (see 
[BGP], |KL] . |KS] . |KuP] ....). In [B G P] the magnetic Schrodinger operator on the planar 
square lattice was considered. The spectrum of operator is expressed in terms of the Lya- 
punov function of the corresponding Hill operator. In [P] these results are extended to the 
case of an arbitrary connected 3D graph with identical edges and the same even potential. 

Schrodinger operators on nanotubes have attracted a lot of attention recently (see |KL] , 
|KET] . [KuP] ....). Authors of [KuP] consider the case of the zigzag, armchair and chiral 
nanotubes with even potentials q. They show that the spectrum of the operator, as a set, 
coincides with the spectrum of the Hill operator. Note that these results can be obtained from 
the results of [P] (K.Pankrashkin, non-published). In some sence the case of even potentials 
q is closed to the case q = 0. Indeed, if q is even, then F_ = [MW] and by Theorem Ol the 
corresponding Lyapunov functions are expressed in terms of the Lyapunov function F for 
the Hill operator. Using the identity F(X) = cosk(A), where k(A) is a quasimomentum for 
the Hill operator, the Lyapunov function is expressed in terms of cosk(A), which is similar 
to the case q = 0, where F(X) = cos y/\. 

The detailed descripiton of spectrum of Schrodinger operator and magnetic Schrodinger 
operator on zigzag nanotube with arbitrary potentials was given in the articles [KL| . [EXT] . 
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|K1] . The operator was represented as the direct sum of quasi one-dimensional operators. 
Research of the spectrum of these quasi one-dimensional operators is based on the analysis 
of the monodromy matrix (see |BBKj . [CKj ). Authors of |KLj describe all eigenfunctions 
with the same eigenvalue. They define a Lyapunov function, which is analytic on some 
Riemann surface. On each sheet, the Lyapunov function has the same properties as in the 
scalar case, but it has branch points (resonances). They prove that all resonances are real 
and they determine the asymptotics of the periodic and anti-periodic spectrum and of the 
resonances at high energy. They show that there exist two types of gaps: i) stable gaps, 
where the endpoints are periodic and anti-periodic eigenvalues, ii) unstable (resonance) gaps, 
where the endpoints are resonances (i.e., real branch points of the Lyapunov function). They 
describe all finite gap potentials. They show that the mapping: potential — > all eigenvalues 
is a real analytic isomorphism for some class of potentials. 

In |KLlj authors describe how the spectrum depends on the magnetic field. They observe 
so-called localization of the spectrum for some values of the magnetic field, when the spec- 
trum of the operator become pure point. This phenominon is typical for quantum graphs 
and was discussed in many papers (see ref. in [PlJ). 

Korotyaev |Klj considers the effective masses for zigzag nanotubes in magnetic fields. Fol- 
lowing |KLj . |KLlj . he consider so called modified Lyapunov function and the corresponding 
quasimomentum. The modified Lyapunov function is entire, then the quasimomentum is an 
analytic function in C+. As for the case of Hill operator, the quasimomentum is a conformal 
mapping of C + onto the domain Re k > 0, Im k > with vertical slits. Identities and a priori 
estimates of gap lengths, heights of slits in terms of effective masses are obtained. These 
values as functions of magnetic field are described. 

Schrodinger operator on armchair nanotube with arbitrary potentials was considered in 
[BBKL], |BBK1] . In |BBKL] all eigenfunctions of the operator with the same eigenvalue 
are described. The Lyapunov function is analytic on some Riemann surface and has branch 
points (resonances). Some example is considered in this paper, which shows that there are 
real and non-real resonances. The detailed analysis of the spectrum is a subject of the 
article [BBK1] . There exist different types of gaps: i) periodic ( ant iper iodic) gaps, where 
the endpoints are periodic (antiperiodic) eigenvalues, ii) resonance gaps, where the endpoints 
are resonances, iii) r-mix gaps, where one endpoint is an antiperiodic eigenvalue and other 
endpoint is a resonance, iv) p-mix gaps, where one endpoint is an antiperiodic eigenvalue 
and other endpoint is a periodic eigenvalue. Asymptotics of the gaps at high energy are 
obtained. 

We present the plan of the paper. In Section 2 we prove Theorems 11.11 11.21 and Theorem 
!2.2l about the eigenfunctions. In Section 3 we prove the main Theorems 1 1 . 3tTL6l Furthermore, 
in this Section we prove Theorems 13.41 13.61 where we give the detail analysis of the spectrum 
including the multiplicity, description of endpoints of gaps etc. Moreover, in this Section 
we prove an existence of gaps independent on the magnetic field for some specific potentials 
(see Proposition 13.51) . 
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2 Spectrum of the operators H% 

We need the following result. 

Lemma 2.1. Each function o w (t) = {^{v^ + te u/ ),e u/ ) = aj, u = (n,j,k) G Z, t G [0,1], 

W? 2 . 6(^2 - gi) W? 2 . W? 2 

ai = a 3 = — a 5 = a6 = — — sm(p — a) = , 02 = 04 = — — sm la = ——. (2.1) 

_ t: Zi 4t 

Proof. Identity ^(r) = |[e 3 , r], e 3 = (0, 0, 1), r G K 3 yields 

M*) = ^([ e 3,r^ + teu^eu) = ^([e 3 ,rj,e w ) = a w (0) = a w , any i £ [0,1]. (2.2) 

Moreover, the vector field .2/ (r) is parallel to the plane x 3 = (this plane is orthogonal to 
the axis of the cylinder C) and £/(r) is tangential to the cylinder surface x\ + x\ = R 2 . 
We deduce that each projection a w ,w G Z of ^(r) to the edge satisfies and 

d\ = Ct 3 = — O5 = Og, CI2 = 04- 

We have e 0>j;k = r 0J+ljk - r 0Jik ,j = 1, 2, which yields 

a o^=2 ([e3 ' ro ^ ] ' r ^ +1 ' fc) - 
If j = 1, then r ,i,fe = R{c 2k , s 2k , 0), r , 2 ,fc = he 3 + R(cos( k , sin( k , 0) and thus 

bR2 /cosC, sinC, 0\ ^ 6jR2 
ai = ao,i, k = — det I I \ = — sm(( k - (p 2k ) = — sm(fj - a) . 

c 2 fc s 2k 0/ 

If j = 2, then r , 3 ,fc = fre 3 + i?(c 2 fc+i, s 2 fe+i, 0) and thus 

bR 2 (C 2k+1 S 2k+1 0\ ^ 2 6i? 2 

«2 = ao,2,fc = — det 1 = — sm(0 2 fe+i - Cfc) = — sin 2a, 

\cosCfc sin C fc 0/ 

since 02fc+i — Cfc = 77 — /3 + a = 2a, which yields (12. ip . ■ 
Now we will prove Theorem 11.11 and the identities 

TrTWo = 2(9F 2 - F 2 - 1), Tr M k = Tr .M - 4s£, (2.3) 

Tr.M 2 = 72F 2 + ^(Tr_M ) 2 -4, Tr A-l* = Tr.M 2 - 8s 2 . Tr M - 16s 2 k c 2 k . (2.4) 

Proof of Theorem 11.11 and identities (12.31) . (12.41) . (i) Introduce the unitary operator % 
in L 2 (T N ) by (<%ff) u = e ttaiJ f u , f = (f w )wez, where is given by (12.11) . We define a modified 
operator H a = ty*Jf B ty in L 2 (T N ). Then H a is given by (# a i% = -f£+qf m u e Z, where 
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/ G ®(H a ). The domain ®(H a ) consists of the functions / = (f u ) ueZ , UZ)^ez € L 2 (T N ) 
and / satisfies the Modified Kirchhoff Boundary Conditions: 

e ia "f Ul {l) = UM = / WB (0), e ia2 ^ 2 (l) = U 3 (0) = U 6 {0), 
e^U A (l) = UM = e Ml /n-i,6, fc (l), e^U 4 (l) = / n ,i, fc+ i(0) = e -™%-i,5, fe +i(l), ( 2 - 5 ) 

/: 2 (o) - e Mi /:,(!) + 4(0) = o, /: 3 (o) - e ^(i) + /: 6 (o) = o, 
-e iai /: 3 (i)+/: 4 (o)-e wi /;-i, 6 ,,(i) = o, -e M v: 4 (i)+/;.i, fe+ i(o)-e- iai /;-i, 5 , fc+ i(i) = o, 

(2.6) 

for all u>j = (n,j, k),j = N 6 , (n, fc) G Z x Zjv- 

Define the operator <S in by Su = (un, Ui, . . . , mat_i) t , u = (tt n )f G C N . The unitary- 
operator 5 has the form S = s k Vk, where Sek = s k ek, tk = -t-(1> s_fc ) s~ 2fc , s ~ kN+k ) 

is an eigenvector (recall s = e*"^); V^u = ek(u,ek) is a projector. The function / in the 
Kirchhoff boundary conditions (12. 5ft is a vector function / = (f w ),uj = (n,j,k) G Z. We 
define a new vector- valued function f n j = (f n ,j,k)k=i e C^, (n,j) G 2i = Z x N 6 , which 
satisfies the equation — • + g/ nj - = A/ nj -, and the conditions 

e iai /n,l(l) = /n, 2 (0) = / n>5 (0), e ia2 /n, 2 (l) = /n, 3 (0) = / n , 6 (0), 

e Ml /n, 3 (l) = / M (0) = e Ml /n-i, 6 (l), e ia2 5/ M (l) = /„,i(0) = e-^/^l), (2.7) 

e jai /:,i(i) - /; 2 (o) - /; 5 (o) = o, e -v; 2 (i) - /; 3 (o) - /; 6 (o) = o, 

e iai /; 3 (l) - fU°) + emi /n-i, 6 (l) = 0, e^Sf nA {l) - /^(O) + e--V;-i, 5 (l) = 0, (2.8) 

for all n G Z, which follow from the Kirchhoff conditions (CD). The operators S and J#b 
commute, then we deduce that J^bPu is unitarily equivalent to the operator H% acting in 
L\T l ) and H a k is given by (H a J) a = -/£ + q{t)f a , where (/ Q ) a6 ^, (/*) oe ^ G L^T 1 ) and 
components f a satisfy the boundary conditions (11.41) . (I1.5p . Thus M'b is unitarily equivalent 
to the operator H a = ©f#£. 

Proof of (ii) and ( 12.31) . ( 12.41) repeats the arguments from [BBKLJ. ■ 

Define the subspace H%(\) = {ip G : H^if) = \ip] for (A, k) G <7 W (#£) x Z N . 

Theorem 2.2. Lei (a, A, fc) G R x a D x Z^. T/ien 

i) Every eigenfunction of7i%(X) vanishes at all vertices ofT 1 . 

ii) There exist the functions ip^ ^ = (tpa ) a eZ! G H%(\), v G N 2 on T 1 such that supp^ '^ 

C U<=o,i(U j6 N 6 r^), each junction = (V»mj )(m^)e2i = {'^m-lj)(mj)ez 1 e Wg(A),n G 

Z. Moreover, each f G ^(A) satisfies: 

f= Yl fnA {n,V \ (fn,ljn,2)n & ef®£ 2 , (2.9) 

(n,i/)eZxN 2 
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Figure 4: The supports of eigenfunctions: a) tp^ ' 1 ^, b) -0( ' 2 ) for (s^, 2 ) 7^ (0, 1), c) ^^ 0,1 \ d) 
.0(o>2) f or 2 ) = (0, 1); ipn'P = CnjYi where each C^- is written near the corresponding 
edge 



where if (sk, 4> 2 ) — (0, 1), 4> = then 



*f(sk,<p 2 ) ^ (0,1), then 



7 z; 5 (o) + s k ^^r n M ? 

Jn,l — , Jn,2 



fn,l = fn,l(0), /n, 2 = /; 2 (0), 

e^f' 6 (0) + ^ a f'(0) 



x x x 2 



x x x 2 



(2.10) 



(2-11) 



K\ = 1 — s k e 2ia (j) 2 7^ 0, x 2 = 1 — s fc e 2ja 4 7^ 0. Moreover, the mapping f — > (f n ,i, fnS)n& is 
a linear isomorphism between Til (A) and £ 2 © I 2 . 



Proof. Proof of i) repeats the arguments from |KLj . 

ii) We will define the eigenfunctions of H% (see FigHJ). If (sk,4> 2 ) 7^ (0,1), then x x 7^ 0, 
x 2 7^ and the function 0)(°> 1 ) is given by 



^ 1) = 0forall(n,j)G(Z\{0,l})xN 6 , and 0^ = 0, j = 1,3,4, ^f> = 0, j = 5, 6 



(0,1) 



■,(0,1) 



0^) ^) 



~ / (0,1) / iu 



-iai 



(0,1) _ ^0,6 



(2.12) 



the function -0(°' 2 ) is given by 

<f = 0, forall(n,j)G(Z\{0,l})xN 6 , and $f = 0, j ^ 4, 4° 2 2) = S fc 0(0 2 - l)e*V, 



^ (0 ' 2) ^ (0,2) ^ (0 ' 2) 
ro,4 _ ^0,3 yo,i 

"~ e~ iai ~ s k (b 2 e ia2 



W,(0,2) ,(0,2) ,(0,2) 

ri,4 ^0,6 _ ^0,5 

e _Wl s k d>e m 



x 2 <^. 



(2.13) 
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If (sfe, <p 2 ) = (0, 1), then x x = x 2 = and the functions ^°' u ^ are given by 

1>£? = 0, (n,j) G (Z \ {0, 1}) x N 6 , VS" 5 = 0, J G N 6 \ {4}, 

<i 1} = = -e- ia 4T = e~ ia 4T = -4T = = ^ (a+ol V, 4T= o, 

<i 2) = <s 2) = <4 2) = 0, < 2 2) = -< 5 2) = e->!° 4 2) = e-0<f = <p. (2.14) 

Using we deduce that tp (0 ' v) satisfy the Kirchhoff conditions COD, <|TT5j> . Thus ^ (0,!y) 

are eigenfunctions of The operator H% is periodic, then each ip( n ' u \ (n, 6 Z x I 2 is 
an eigenfunction. We will show that the sequence ip( n ' u \ (n, v) G Z x N 2 forms a basis for 
Wg(A). Let h = E(n,, )e zxN 2 «„,^ (ni!y) = 0. Then ©, © imply for 2 ) ^ (0, 1) 

= ^Ir„ l8 = («n,i + a ni2 s fc 0e ia )x 2 ^, = /i|r n , 6 = -(a„ 5 i0e ia + a n)2 )x 2 e- ia V, n G Z, 

which yields a nji , = for all (n, v) G Z x N 2 . Hence ?/;( ri ' I/ ) are linearly independent. The 
similar arguments show that ift( n > u > are linearly independent for (sfc,0 2 ) = (0, 1). 
For any / G T~Ct(\) we will show the identity ( 12.91) . i.e., 

/ = /, where / = ^ /n,^ (n '" } and f n , u are given by (J2H). (2.15) 

(n,i/)eZxN 2 

/From A G cx D , we deduce that f\ VnJ = fn d (0)<p. If (s fc , <p 2 ) ^ (0, 1), then identities ([2?TT|) -(|2]) 
provide for all n G Z 

/k, 6 = £ Arf^k. = (/n,i + l,2S k <Pen^ = f n ,M* = /k«, 
(£,i/)eZxN 2 

/k >6 = S,^ (V V„, 6 = -(7nA^ a + 7n, 2 )x 2 e- ia V = /;6(0)^ = /|r„, 6 . (2.16) 
If (s fc , 2 ) = (0, 1), then identities (12101) .(121) give for all n G Z 

/| rnil = = fU )f = /kl, 7r„, 2 = /n,2^ = /^(O^ = /|r n , 2 - (2.17) 

Identities (ESD, (I2TT1) yield £ |/„^| 2 < oo and / G L 2 ^ 1 ), since / G L 2 ^ 1 ). 

Note that / satisfies the Kirchhoff conditions (11 .41) . (11.51) and — f£ + qf a = X f a , a G Z\. 
Consider the function u = f — /. The function u = at all vertices of T 1 and then 
«nj = C nJ ip, (n,j) G Z x N 6 . If (s fc ,0 2 ) 7^ (0,1), then (I2.16j) gives C„ i5 = C„ >6 = 0. Let 
n G Z and assume that C n> i = C. Then the Kirchhoff boundary conditions (11.4l) - (11.5j) 
yield C n>2 = C0e iai , C n , 3 = C0V a , C n>4 = C0V( a+ai ), C n ,i = Cs k ^e 2ia . Since C n>1 = C, 
we obtain C = and u = 0, which yields (12151) . If (s fc ^ 2 ) = (0,1), then (12171) gives 
Cn,i = Cn,2 = 0. The Kirchhoff boundary conditions (ll.4p - (ll.5l) yield C n $ = 0, and then 
C n ,A = 0,n G Z. Assume that C n)3 = C. Then C ni6 = — C 0i 3 = — C and C n+ i j3 = — C„ i6 = C 
for all fiGZ. Due to w G L 2 (r x ), we have C = and u = 0, which yields ( 12. 151) . 
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The mapping / — > (/ n ,w)(n,j/)eZxNa i s a linear and one-to-one mapping from W£ onto 
I 2 © I 2 . Then it is a linear isomorphism. ■ 

Proof of Theorem 11.21 (i) Using the arguments from [BBK],[KL] and identities (11.101) . 
(ESD, (E3D we obtain (TTT2]) . (Q51) . 

(ii), (iii) Proof of the statements and identities (11.151) repeats the corresponding argu- 
ments from [KEj. Identities (021) . fTTTTD yield —a) = F N - ktV (-,a), F ki „(; a + f ) = 
a), , A;) G N 2 x Z w . Then identities Q5| yield (OS1) . By Theorem^ each point 
Ao G a pp(Hk) nas an infinite multiplicity. ■ 

Lemma 2.3. Le£ (a, fe) G K x Zyv- TTien i/iere exists an integer n > 1 stzc/i that 

(i) The function = Dq = det(AAk + h) has exactly 4n zeros, counted with multiplicities, 
in the domain {A : \VX\ < ttuq} and for each n > uq, exactly two zeros, counted with 
multiplicities, in each domain {A : |vA — 7m — | ± arcsin || < |}. There are no other zeros. 

(ii) Each function D\ = det(.Mfc — h),k G Z N has exactly 4n + 2 zeros, counted with 
multiplicities, in the domain {A : \VX\ < ttuq + f }, and for each n > no exactly one simple 
zero in each domain {X : \y/~X — nn ± arccos ^- 5 ~ 4cfc | < |} ; exactly one simple zero in each 
domain {X : |\/A — nn ± arccos ^ 5 g 4cfc | < y }, if Sk ^ 0, and exactly two zeros, counted with 

multiplicities, in each domain {X : |\/A — im ± arccos ^ 5 g 4cfc | < ~} ; if s k = 0. There are no 
other zeros. 

(iii) Let c k ^ 0, Sk ^ 0. Then the function p k has exactly 2n zeros, counted with multiplic- 
ities, in the domain {X : \VX\ < 7m }, and for each n > n exactly one simple real zero in 
each domain {X : \y/\ — (nn — | ± arcsin y)| < y }. There are no other zeros. 

Proof. Proof uses the Rouche theorem and repeats the arguments from [KL]. ■ 
Below we need some identities for Df. Substituting (jl.121) into (11.181) we obtain 

D+ = (9F 2 - g h>1 )(9F* - g k , 2 ), D k = D = (9F 2 - h)(9F 2 - h 2 ) on R, (2.18) 

where 

9k,u = 5 + F 2 + (-1Y2^F 2 +Acl h v = (1 + (-ir|F_|) 2 , (u, k)eN 2 x Z N . (2.19) 
Lemma 2.4. Let (a, k) G R x Z N . Then 

(i) For all [v, n) G N 2 x N the periodic and antiperiodic eigenvalues satisfy 

A'liW^^tH' X k u f(a) = X°^(0), 7 '»=7° P (0), T ° p = x n , p = 2n - 1, 

(2.20) 

< »M < A i£ < niin{Ai';~ +1 , A°J } ^ max{A^ +1 , A°;J} ^ A 2 ; p+1 ^ A", (2.21) 

A° 2 ;; < Aj;+ i ^ u k (Xl p ) > ; A°£ > A^+i & u k (X° 2 +) > 0, (2.22) 
Un>i7° 2 n-i = {A G R : 9F 2 (A) < fc„(A)}, U n >ol k u ,2n = {A G R : 9F 2 (A) > ^(A)}, (2.23) 
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Figure 5: Functions 9F 2 ,g k t ,h ± and fk 



where 

7 Jo = (-00, Aj+), 7 J„ = (A*" A*;+), (u, n)EN 2 x N. 
^ If s\(a) < sf(ai) /or some (ai,£) 6lx Z^v, £/ien /or all n ^ 1 

A2,'2n-2( a ) < A2,2n-2( a l)) -V,2n-2( a l) < ^l,'ii-2( a ) > ^1.2n ( a ) < ^1.2n ( a l) 7 A 2 ' i2 n( a l) < ^2,2n( a ) ' 

(2.24) 

Proof. Proof uses identities (I2.18P and properties of the function F and repeats the argu- 
ments from |BBKlj . Lemma 2.2 (see also FigJS]). ■ 

3 Proof of Theorems 

Let Rk = {A G R : pfc(A) > 0}, if p k ^ 0, and Rk = R, if = 0. In Lemmas I3.1H3.3I we 
describe the set o^ v = {A G R : F ku {\) G [—1, 1]} in terms of F. 

Lemma 3.1. Let a G R. For all k G Zjv arid A G i/ie following relations hold true: 

F ktV (X)<l iff 9F 2 (A) < o^(A), ^ = 1,2, (3.1) 
F M (A) > -1 ^ {9F 2 (A) > MA) ^ |F_(A)| < c 2 }, (3.2) 
F k>2 (\) > -1 ^ |9F 2 (A) > h 2 (X) or {9F 2 (A) < ^(A) and |F_(A)| < c 2 } j. (3.3) 
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Proof. Proof repeats the arguments from |BBK1] . Lemma 3.1. Using the identities (11.121) . 
(I2.19P we write the functions F k>v in terms of g k)U , h u and prove (I3.1I) - (I3.3I) . ■ 

Now we describe the zeros of p k and the functions F kjl/ on the interval x n . 

Lemma 3.2. Let ^ for some (a, k) G R x TL^. Then 

(i) The following relation holds true: 

R \ R k C U„^ix re , where x„ = (A^n-i, ^i,t n -i)- ( 3 - 4 ) 

All real zeros of pk belong to the set U n ^ix n . Each interval x„,n ^ 1 contains even number 
^ of zeros of pk, counted with multiplicities. 

(ii) Let x n <f_ Rk- Then for each a = ± the following relations hold true: 

signt>fc = const on each x kn C Rk, (3.5) 

if Vk(X) > for some A G x£ n , t/ien F k ^ < F k ,\ < —1 on x kn , (3.6) 

u fc(A) < for some A G x£ n , t/ien — 1 < F fe2 < F k> i on x£ n , (3.7) 

where 

*n,k = ( X l,2n-n r k,n)> K,k = ( r Jn> ^l^-l). n > 1 ( 3 - 8 ) 

// x n C Rk, then 

Vk > and Ffc i2 < -Pfc,i < —1 on x n . (3-9) 
(m,) If Vk{X) < /or some A G x^ ; t/ien has even number ^ 2 of zeros on TQ. 

Proof. Proof repeats the arguments from [BBKlj, Lemma 3.2. The last identity in (I1.12p 
gives 

p = cl(9F 2 -fk), where f k = s 2 k (l - ^ j. (3.10) 

Identities (12.191) yield hi ^ f k . Using the properties of F we obtain (13 .4p (see also FigEJ) • 
Using (13.101) . Lemma 13.11 and the more detail analysis (see [BBKlJ) we obtain (I3.5l) -( l3~9l 
and the statement (hi). ■ 

For each (k, v) G x N2 we introduce the sets 

©V = U A S~-i> X %U X °4 , A";;+i]) , 6f = |J xf^, AT ± = G N : u fc (A°£) < 0}, 

n^l ^ ' o-=±,neiV CT 

(3.11) 

where p = 2n — 1. The set & kjl/ is a "stable" part of cr k>v = {A G R : F k , v (X) G [—1, 1]}, and 
6f is an "unstable" part of a Kv (see (I3T3D . fl3TT4l) ). 

Lemma 3.3. Let {a,k,v) G R x Z w x N 2 . T/ien i/je following identities hold true: 

G k , v = {A G R : h v {X) < 9F 2 (A) ^ <MA)}, (3.12) 

& ^ = \{XeR k :9F 2 (X)^h 1 (X) andv k (X)^0}, if c k ^ 
'"• '/>/, 

o-fc^ = &k, v U 6f • (3.14) 
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Proof. Proof uses the results of Lemma 13.11 aud repeats the arguments from |BBK1] . Lemma 
3.3. Relations (jQHD and Lemma O yield fl3~T2]) . (l3~T3j) . Then Lemma O gives <KW\ . m 
Proof of Theorem 11.31 Substituting (13. lip into (13. 14j) and using the identity a ac (H%) = 
CM U a K2 we obtain (fl~20l) . (OB . Identities A*;J(a) = A°;J(0) are proved in (l2~20D . Using 
(fT211 and the definition of x n in (fl~23l) we obtain (Q2l) . ■ 

Theorem 3.4. (multiplicity of the spectrum) Let (a, k, n) E [0, ^] x Z^r x N,p = 2n - 1. 

T/ien t/ie following relations hold true: 

E k 4 = E° 2 £, ^ min{E k 2 >-, E^} < max{E k 2 >~ , Ej£.J ^ E k '~ 

^ E*'+ < min{*&i> 4'+} < max{4'- +1 , £ 2 fc p +} < £ 2 fc ;; +1 , (3.15) 

4" < «*(^-)>0; 4; + i< «*(^i + )>0. (3.16) 

Moreover, 

(%) If 7^ 0, t/ten t/ie spectrum of H% in &^ has multiplicity 4- 

(ii) If &k = &k,i H 6fc,2 7^ #ien t/ie spectrum of H% in &t has multiplicity 4- 

(Hi) The spectrum of H% in a ac (H%) \ (&k U (5f) has multiplicity 2. 

Proof. Relations (OB and (ELD) yield (13151) . Relations (12T22D imply (l3~T6l) . Proof of the 
other statements repeats the arguments from [BBKlJ, proof of Theorem 1.2. Relations (13. 7p 
provide the statement (i), identity (I3.14p yields the statements (ii), (iii) (see [BBKlj ). ■ 
Remark. 1) Using identities (fL"2T|) . fl3~TD) we deduce that 6£ ^ iff {E^ = r% >n , x£ n ^ 

for some (n, a) G N x {+, — }}. 

2) Identities (fl~2"ll and (131TD show that 6 fc ^ iff > E k '^ x (or > for 

some n ^ 1. Relations (13.161) yields Uk(E 2 '~) > (or Uk{E 2 p) > 0) for this n. Then the 
spectrum of if£ in the interval {E x '^_^ E 2 ~) (or (E 2 p, E k ] p+1 )) has multiplicity 4. 
Below we need the asymptotics from [M] 



F(A) = cosv / A+ 1 J , F_(A) = ^ e j as |A| ^ oo. (3.17) 

v|A| v|A| 



Proof of Theorem 11.41 (i) Estimates (I3.15P show that the intervals G u ^fi = ( — °°^E k, Q). 
G Vlk , n = E%), {y, n) e N 2 x N, satisfy: 

G lt k,o n G 2 ,fc, m = for m £ {0, 1}, Gi,fc,2n-i n G^m = for m ^ 2n - 1, 

Gi^^n n G 2 ,fc,m = for m <=£ {2n - 1, 2n, 2n + 1}. 
Then the gaps n , n ^ in the spectrum if£ satisfy 

^fc,0 = ^l,fe,0 H G2,k,0, G\ 2n = Gl,A:,n H G^^k^ni 

^fc4n-3 = Cl,fc,2n-2 H G^fc^n-lj Cfc 4n-l = ^l,fc,2n H G 2 ,fc,2n-1 > (3.18) 
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n ^ 1, which yields all identities in (ll.24p . Estimates (I2.2ip give all inclusions in (11.241) . 
Relations fl2~2^j) give (Q5]) . 

For (a, fc) = (0, 0) the first identity in (OBj) is proved in [BBKlj . If a ^ 0, then s fc ^ 0. 
Let Cfc 7^ 0. Asymptotics (13. 17f) implies -F-(A) — > as A — > +oo. The last identities in (11.211) 
give i?-^ = r^ n , n ^ no- By Lemma [2731 (iii), r A T n < r£ n for such ra. Hence the intervals 
Gl 4n _ 2 7^ in this case. Let Ck = 0. The last identities in (11.211) give E^ p = A^p,n ^ 1. 
Identities f l3~TH . fl3~T2l yield A { are zeros of 9F 2 — hi, where hi is given by (12.190. We 
have hi(X) — > 1 as A — ► +oo. The properties of F show that the function 9F 2 — hi has 
only simple zeros on (Ao, +oo) for sufficiently large Ao > 0. Then A^~ < X{p for all n ^ n . 
Hence G^ 4n _ 2 ^ in this case, which yields the first identities in (11.261) for (a, k) ^ (0, 0). 

If (a, k) 7^ (0, 0), then s k ^ and estimates (13.161) give E^ An _ 3 > E^ An _ 3 and -E A T 4n _ 1 > 
E£ 4n _ v n > n . Hence 2n _ 1 = for such k, n, which gives the second identities in (11.261) . 
(ii) If ffc(A) > on x n , then identities (ll.2ip show that E^(a) = A^ p . Then (11.241) gives 

Proof of Theorem II 51 Using the results of Theorem 11.21 (iii) we obtain (ll.27p . 

(i) Theorem O (i) provides a ac (H a ) = R\ U n>0 G«, where gap G a n = n fceZjv G^ n . The 
first relations in (11.251) show G\ n = GQ An ,n ^ 0. The second relations in (11.251) imply 
G% n -\ = G a mi 2n -\-. n ^ 1; where mi = | or The relations G a kAn _ 2 C x n give G% n _ 2 C 
K n . The relations 7„_i C Gl 4n ,r] n G [E 1 ' 2n _ x > E x 2n-i] gi ve the corresponding relations 
7n C G4 n ,?7 n G [-E~2„-i) ^i~2n-i]- Thus, we have proved all relations in fll.28j) . Relations 
fOHjl yield (TT7291 . Relations (Ofil) yield (OOl) . 

(ii) In order to prove asymptotics (II. 311) we assume that J* q{t)dt = 0. Let m — 0. Using 
the first identities in (OTjl and ffl~24l) we obtain (E An , E+J = G% n = G a QAn = {E 2 ' 2n , ) = 
(A^2n5 A°^ n ). Identities (12.191) and Lemma [2731 (ii) show that Ef n = X 2 ^ n are zeros of the 

function 9F 2 — g ,2 and a/ ~E~f n = nn ± 9 + , where |e^| ^ | for large n. Let A = Ef n and 
e = e±. Then ( 13717]) give 

9F 2 (A) = 9 cos 2 (vm ± # + e) + 0(n~ 2 ) = 5 + 4|c | =F 9e sin(2# ) + 0(e 2 ) + 0(rr 2 ), 

and F_(A) = 0(n _1 ) as n -> +00. Identities f |27l9j) imply # ,2(A) = 5 + 4|c | + 0(n~ 2 ). 
Then 9F 2 (A) - #0,2 (A) = T9£sin(2# ) + 0(e 2 ) + 0(n~ 2 ). Hence e = 0(n" 2 ) and ^/E~f n = 
irn ±6 + 0(n~ 2 ), which yields (ll.3ip for m = 0. 

Let m = 1. Identity G£ = n fceZjv G^ n and p.24p impliy E 4n _ 2 = max teZjv -Ef^i, 
Etn-2 = m ^kez N Ei2 n _i. Identities (11.121) show that for large n ^ 1 the energies E^ 2 i: n _ 1 = 
r fc2n are zeros °f the function 9F 2 — where fk = s\(l t)- By Lemma [2~31 (i), (iii), 

' k 

\j Ei'^n- 1 — 71-72 — f =t ^fc + ; ®k = arcsin ^, where |e*| ^ ^ for large n. Let A = E^ 2n _ 1 
and £ = e^. Then (I3.17P give 

9F 2 (A) = 9 cos 2 (t™ - I ± 9 k + e) + 0(n" 2 ) = s 2 ± 9e sin(20 fc ) + 0{e 2 ) + 0(r?r 2 ), 
and F_(A) = ^^(n- 1 ). Then / fe (A) = s 2 + 0(n~ 2 ) and 9F 2 (A) - / fe (A) = ±9esin(2^) + 
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0(e 2 ) + 0{n- 2 ). Then e = 0(n~ 2 ), which yields E^_ x = (vrn - § ± 9 k + 0(n- 2 )) 2 . This 

asymptotics and the identity 9\ = minkez N 9k give ( jl.3ip for m = 1. 
(iii) Theorem 11.41 (ii) yields the statement. ■ 

Proof of Theorem [Tm (i) If q e L 2 even (0, 1), then F_ = 0. Let c k = 0. The last identity in 
(I1.12p implies pk = and identity (11. 19ft yields Vk = 0. Using identities (11.211) = A$'^ l _ 1 
and (Q3l) . (Q4l) provide G^ 4n _ 2 = x n . Let c fc ^ 0. Then (091) yields w fe = -c 2 , < 0. 
Identities (11.121) show that the resonances r^ n are zeros of 9F 2 — s\, hence they are real. 

Identities (ESQ) imply = r% n and JOi 'gives G% M _ 2 = (r~ n , r+J. 

The identities (031 for ( a > fc ) = (0,0) are proved in [BBKlj . Let (a, k) ^ (0,0). Then 
Sfc(a) > and properties of the function F yield r^ n < r£ n . Hence G^. An _ 2 ^ 0, which yield 
the first identity in (11.331) for this case. Moreover, (II . 191) shows Uk ^ and estimates (13 . 161) 
S ive E kAn-s > and ^4n-i > ^4n-P 11 ^ E Hence G^ 2n _ 1 = 0, which gives the 

second identities in (ll.33p . If s k (a) < Sf(ai), then the properties of the function F yield 
±r*>i) < ±r|». Then (TOl gives G?, te ^„., 

(ii) The first identity in (11.351) was proved in [BBKLJ. The other relations (ll.35j) and (II .34j) 
follow from the statement (i). We prove (PHI) . (PT1) . Recall that E± n (a) = A 2 ' 2n (a) are 
zeros of the function 9F 2 — g ,2- Since F_ = 0, identities ( 12.191) imply g ,2 = 5 + 4 cos a. Then 
9F 2 — ^0,2 = / + 4(1 — cos a), where / = 9F 2 — 9, and Ef n (0) = A^ (see the first identity in 
(I1.35P ) are zeros of the function /. Let 7 n ^ 0. Then 

/(A±) = 0, f(\i) = 18F(A±)F'(A±) = -M± ^ 0, n > 1, 

where we have used the identities = — F(A^)F'(A^) (see [KKJ). Let Ef n (a) = X^+e,e = 
e ± . Then f(Ef n (a)) = e/'(A±) + 0(e 2 ) = -sM± + 0(e 2 ), and using the asymptotics 
4(1 — cos a) = 2a 2 + 0(a 4 ),a — > 0, we have 

= 9F 2 (£^ n (a)) - g 0>2 (Et n (a)) = -18M±e + 2a 2 + 0(e 2 ) + 0(a 4 ), a -> 0, (3.19) 

which yields e = 0{a 2 ). Substituting this asymptotics into (I3.19P we obtain e = -f^+0(a 4 ), 
which yields ( 11.361) . 

Let 7 n = 0, i.e. E = A~ = A+. Then 

/(E) = /'(£) = 0, /"(E) = 18F(E)F"(E) = -18\F"(E)\ < 0. 

Let E± n (a) = E + e,e = e ± . We have f(Ef n (a)) = ^f'(E) + 0(e 3 ) = -§\F"(E)\s 2 + 0(e 3 ) 
and 

= 9F 2 {Ef n (a)) - ^ 0)2 (^ n (a)) = -9|F"(£)|e 2 + 2a 2 + 0(e 3 ) + 0(a 4 ), (3.20) 

which yields e = O(a). Substituting this asymptotics into (I3.20p . we obtain e 2 = 9 ^ 2 f^ + 

0(a 3 ), which yields ( Q71) . 

We prove (Q8|) . Using G£ 4n _ 2 C G^ 4n _ 2 and (OH . (ESP we obtain Ef n _ 2 {a) = 
-^i,2n~i( a ) = r o,n( a )- Identities (11.12)) show that r^ n (a) are zeros of the function 9F 2 —sin 2 a 
and r o.n(°) = r o,n(°) = Vn- Let r± n (a) = Vn + £ ,e = e± Then = 9F 2 (rJ n (a)) - sin 2 a = 
9(F'(7] n )) 2 e 2 -a 2 + 0(e 3 ) + 0(a 3 ), which yields £ = O(a). Then e = ±^^ + 0(a 2 ), which 
gives (OKI) . ■ 
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Proposition 3.5. Let q = q £ = -5(t — § — 2e),e > 0. Then for any n G N there exists 
Eq > suc/i i/ia£ 

Gl n . 2 {q £ ) = G° An _ 2 (q £ ) ± a// (a )£ ,n) 6lx (0,£„) xN„ r (3.21) 

Proof. We have F_ (A, g e ) = S1 ° £ 4 ^^ = 2 + 0(£ 2 ), e — > uniformly on any bounded subset of 

C (see, for example, [BBKL] ). Then |F_(A,g e )| ^ cos 2 a for each (a, e, A) G Ex (0, e ) xC (R) 
for any i? > and some £o > 0, where Co(R) = {A G C : |A| < i?}. Using Theorem 11.51 (hi) 
we obtain G An _ 2 (q £ ) = G° An _ 2 (q £ ) = x n (q £ ) for all (a,e,n) G M x (0,e ) x N no for any n GN 
and some e > 0. Then identities Q| yield (l3~2Tj) . ■ 

In order to describe the gaps in the spectrum of H a we need 

Definition 1. Let g = (Ai, A2) be a gap in the spectrum of H% or H a . 

(i) If Ai, A2 are zeros of {or D k ), then g is a periodic (or antiperiodic) gap. 

(ii) If Xi, A2 are zeros of pk, then g is a resonance gap. 

(Hi) If one of the numbers Ai, A2 is a zero of and other is a zero of (or pk), then g 
is a p-mix gap (or r-mix gap). 

In our armchair model there is no a gap (Ai, A2), where one of the numbers Ai, A2 is a zero 
of D£ and other is a zero of pk- 

Theorem 3.6. Let (a, k,n) G [0, x Z^r X N. Then (i) G\ An are periodic gaps, 
(ii) G^2n-i are p-mix gaps. Let 7^ 0. Then G% 2n -i = ®> n ^ n o f or some n o ^ 1- 
(Hi) G^ An _ 2 are antiperiodic, or resonance, or r-mix gaps. Moreover, for some n ^ 1 

if s k = 0, then G^ An _ 2 are antiperiodic gaps or G% An _ 2 = 0, and G% An _ 2 = for n ^ n 
if c k 7^ 0, Sk 7^ 0, then G\ An _ 2 are resonance gaps for n ^ n 
if Ck = 0, then all G% An _ 2 are antiperiodic gaps. 

Let, in addition, q G L>l ven (0, 1). Then for all n ^ 1 

if s k = 0, then E^n-i = = Vn and G a kAn __ 2 = 

if c k 7^ 0, Sk 7^ 0, then G\ An _ 2 are resonance gaps. 

Proof. Identities ( ll.2ip . (11.24D show that G k An are periodic gaps, G k 2n _ x are p-mix gaps and 
G kAn _ 2 are antiperiodic, or resonance, or r-mix gaps. Asymptotics f)3. 17j) implies F_(X) — > 
as A — ► +00. Estimates f)3.16p give that E kAn _ 3 > E^ An _ 3 and E^ n _ x > E^ An _ l for k 7^ 
and large n > 1. Hence G\ 2n _ 1 = for such k, n. 

If Sk = 0, then ( 11.12ft gives pk = §F 2 c\ and r~ kn = r kn = r\ n . Identities ( ll.2ip give 

{^o,± ^ (\ 0,± ) > 
l ' P J ,±V ^ • Hence G Mn-2 are antiperiodic gaps or G% M _ 2 = 0. 

rj n if ^(A^ ) < 

Moreover, Vk(X^) < and = E\'^ = r\ n for all large n ^ 1. Hence G^ 4 „_ 2 = for 
sufficiently large n ^ 1. If Ck = 0, then ffc(A) ^ for all A G M and identities ( 11.211) provide 
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E x ' p = . Hence G k \ in _ 2 are antiperiodic gaps. Let c k 7^ 0. Then identities (11.211) give 
G% An _ 2 are antiperiodic, or resonance, or r-mix gaps. Using F_(A) — > as A — >■ +00 we 
obtain ffc(A) < for large A > 0, then G k 4n _ 2 are resonance gaps for sufficiently large n ^ 1. 

If 9 G Lg„ en (0, 1) and 7^ 0, then F_ = and v k < 0. Identities (11.211) show that 
E^f = r^ n in this case. The last identity in ffl~12l) yield p k = (9F 2 - s 2 k )c 2 k . If c k ^0,s k ^ 0, 
then the properties on the function F show that r~ kn < r\ n for all n ^ 1. Then G1 in _ 2 = 
( r fcn' r fcn)' n ^ 1 are resonance gaps. ■ 



References 

[ALM] Amovilli, C; Leys, F.E.; March, N.H. Electronic energy spectrum of two-dimensional solids and a 

chain of C atoms from a quantum network model. J. Math. Chemistry 36 (2004), no. 2, 93-112. 
[ALM1] Amovilli, C; Leys, F.E.; March, N.H. Topology, connectivity, and electronic structure of C and B 

cages and the corresponding nanotubes. J. Chem. Inf. Comput. Sci. 44 (2004), 122-133. 
[BBK] Badanin, A.; Briming, J.; Korotyaev, E. The Lyapunov function for Schrodinger operators with a 

periodic 2x2 matrix potential. J. Funct. Anal. 234 (2006), no. 1, 106-126. 
[BBKL] Badanin, A.; Bruning, J.; Korotyaev, E.; Lobanov, I. Schrodinger operator on armchair nanotube. 

I. Preprint arXiv:math/0707.3909, 2007. 
[BBK1] Badanin, A.; Bruning, J.; Korotyaev, E. Schrodinger operator on armchair nanotube. II. Preprint 

arXiv:math/0707.3909, 2007. 
[BGP] Bruning, J.; Geyler, V.; Pankrashkin, K. Cantor and band spectra for periodic quantum graphs with 

magnetic fields. Comm. Math. Phys. 269 (2007), no. 1, 87-105. 
[Ca] Carlson, R. Spectral theory and spectral gaps for periodic Schrodinger operators on product graphs. 

Special section on quantum graphs. Waves Random Media 14 (2004), no. 1, S29-S45. 
[CK] Chclkak, D.; Korotyaev, E. Spectral estimates for Schrodinger operators with periodic matrix po- 
tentials on the real line. Int. Math. Res. Not. (2006), Art. ID 60314, 41 pp. 
[Co] Coulson, C. A. Note on the applicability of the free-electron network model to metals. Proc. Phys. 

Soc. A 67 (1954), 608-614. 
[Ha] Harris P. Carbon Nanotubes and Related Structures. Cambridge Univ. Press., Cambridge, 1999. 
[K] Korotyaev, E. Conformal spectral theory for the monodromy matrix. Preprint arXiv: math/0612284, 

2006. 

[Kl] Korotyaev, E. Effective masses for zigzag nanotubes in magnetic fields. Preprint 

|arXiv:math-ph/0703017] 2007. 
[KK] Kargaev, P.; Korotyaev, E. The inverse problem for the Hill operator, a direct approach. Invent. 

Math. 129 (1997), no. 3, 567-593. 
[KL] Korotyaev, E.; Lobanov, I. Schrodinger operators on zigzag nanotubes. Ann. H.Poincare, 8 (2007), 

1151-1176. 

[KL1] Korotyaev, E.; Lobanov, I. Zigzag periodic nanotube in magnetic field. Preprint arXiv: 
math/0604007, 2006. 

[KS] Kostrykin, V.; Schrader, R. Quantum wires with magnetic fluxes. Comm. Math. Phys. 237 (2003), 
no. 1-2, 161-179. 

[Ku] Kuchment, P. Graph models for waves in thin structures, Waves in Random Media, 12 (2002), 
R1-R24. 



[KuP] Kuchment, P.; Post, O. On the spectra of carbon nano-structures. Preprint arXiv:math-ph/0612021 
2007. 

[MW] Magnus, W.; Winkler, S. Hill's equation. Dover Publications Inc., New York, 1979. 
[M] Marchenko, V. Sturm-Liouville operator and applications. Basel: Birkhauser, 1986. 
[Mo] Montroll, E. W. Quantum theory on a network. I. A solvable model whose wavefunctions are ele- 
mentary functions. J. Math. Phys. 11 (1970), no. 2, 635-648. 



22 



[P] Pankrashkin, K. Spectra of Schrodinger operators on equilateral quantum graphs. Lett. Math. Phys. 
77 (2006), no. 2, 139-154. 

[PI] Pankrashkin, K. Localization effects in a periodic quantum graph with magnetic field and spin-orbit 

interaction. J. Math. Phys. 47 (2006), no. 11, 17 pp. 
[Pa] Pauling, L. The diamagnetic anisotropy of aromatic molecules, Journal of Chemical Physics, 4 

(1936), 673-677. 

[RS] Rucdcnbcrg, K.; Scherr, C.W. Free-electron network model for conjugated systems. I. Theory, The 

Journal of Chemical Physics, 21 (1953), 1565-1581. 
[SDD] Saito, R.; Dresselhaus, C; Dresselhaus, M. Physical properties of carbon nanotubes, Imperial College 

Press, 1998. 



23 



